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ABSTRACT 

Supernova remnants are thought to be the dominant source of Galactic cosmic rays. This requires 
that at least 5% of the available energy is transferred to cosmic rays, implying a high cosmic-ray 
pressure downstream of supernova remnant shocks. Recently, it has been shown that the downstream 
temperature in some remnants is low compared to the measured shock velocities, implying that addi- 
tional pressure support by accelerated particles is present. 

Here we use a two-fluid thermodynamic approach to derive the relation between post-shock frac- 
tional cosmic-ray pressure and post-shock temperature, assuming no additional heating beyond adi- 
abatic heating in the shock precursor and with all non-adiabatic heating occurring at the subshock. 
The derived relations show that a high fractional cosmic-ray pressure is only possible, if a substantial 
fraction of the incoming energy flux escapes from the system. Recently a shock velocity and a down- 
stream proton temperature were measured for a shock in the supernova remnant RCW 86. We apply 
the two-fluid solutions to these measurements and find that the the downstream fractional cosmic-ray 
pressure is at least 50% with a cosmic-ray energy flux escape of at least 20%. 

In general, in order to have 5% of the supernova energy to go into accelerating cosmic rays, on 
average the post-shock cosmic-ray pressure needs to be 30% for an effective cosmic-ray adiabatic 
index of 7 cr = 4/3. 

Subject headings: acceleration of particles - cosmic rays - shock waves - supernova remnants 



1. INTRODUCTION 

Although there have been many theoretical and obser- 
vational advances in understanding cosmic-ray accelera- 
tion by supernova remnant (SNR) shocks over the last 
decade, it is still not yet clear whether SNRs are capa- 
ble of putting more than 5% of their energy into cosmic 
rays. This number is necessary in order to explain the 
cosmic-ray energy density in the Galaxy, given the size 
of the Galaxy, the supernova rate, and the co smic-ray es- 
cape time (e.g. lGinzburg fc Svrov atskii 19671 for an early 
discussion). 

The theoretical advances consist, among others, of 
a better understanding of magn etic field a mplification 
by cosmic-ray streaming (e.g. I Belli 120041 ). and ma- 



jor improvements in the self-consiste nt modeling of ef- 
ficient cosmic-ra y acceleration (e g. iBlasi et alj 120051 : 
iKang et al]|2009t Iviadimirov et alll2008D . Thesesimula- 
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tions show that SNR shocks can transfer up to ~ 80% of 
the kinetic energy flux entering the shock to cosmic rays. 

On the observational side progress has been made 
through X-ray observations which show that many young 
SNRs are actively accelerating electrons up to 10-100 
TeV, and they show that the magnetic fie l ds must in- 
deed be amplified (e . g.lVink fc Laming! 20031: IVoelk etHI 
120051 : iBamba et all 120051 : iBalletl I2006D . Another ma- 
jor source of progress has been the coming of age of 
TeV gamma-ray astronomy. TeV gamma-ray observa- 
tions have s hown us that many young SNRs are TeV 
sources (e.g. lAharonian et al.1 120011 120041: lAlbert et all 
120071: lAcciari et all 120101 ) . This proves that SNRs are 
capable of accelerating particles up to at least 100 TeV. 

j .vink@astro-uu.nl 



Unfortunately, a correct interpretation of the emission 
mechanism, i.e. pion-decay or inverse Compton scatter- 
ing, is necessary for understanding the fraction of en- 
er gy contained by acce lerated particles (see the review 
bv lffinton fc Hofmand[2009[ ). 

Recent GeV gamma-ray observations with the Fermi- 
LAT even show that, whatever the emission mechanism, 
a bright young SNR like Cassiopeia A (Cas A) has only 
transferred ~ 2% of its explosion energy to accelerated 
particles. The interpretation is not yet clear: are Cas 
A and other young SNRs beyond their peak in accelera- 
tion power, and did most of the highest energy particles 
escape? Or did they not yet reach their full potential 
as sources of cosmic rays? Or are young SNRs not the 
dominant sources of Galactic cosmic rays? 

Recently, another method of measuring cosmic-ray ac- 
celeration efficiency has drawn attention. It consists of 
measuring the proton temperature behind SNR shocks 
that are suspected of being efficient cosmic-ray acceler- 
ators, in view of their high velocity or multi-wavelength 
properties. Without cosmic-ray acceleration there ex- 
ists a simple relationship between the post-shock plasma 
temperature and shock velocity. However, when a shock 
accelerates cosmic rays, less energy is available for plasma 
heating. The proton temperature can be deduced from 
thermal Doppler br oadening of Ha line emission from 
behind fast shocks (jHengl 12010). Note that, unlike the 
electron temperature, the proton temperature is for cos- 
mic abundance close to or equal to the mean plasma 
temperature. 

This method was used by iHelder et al.1 (|2009f ) for an 
X-ray synchrotron emitting shock in the SNR RCW 
86. The X-ray synchrotron emission indicates that the 
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shock is actively accelerating particles to > 10 TeV en- 
ergies. Moreover, RCW 86 is detected by H.E S.S. as a 
TeV gamma-ray so urce (|Aharonian et all 120081 ). Indeed 
iHelder et alJ (|2009( ) found that the plasma temperature 
is a factor of at least three lower than expected given the 
measured shock velocity. However, this value could not 
be directly translated into a fractional cosmic-ray pres- 
sure behind the shock, but was translated into a lower 
limit of fractional cosmic-ray pressure of > 50%. This 
lower limit was derived using the Ra nkine-Hugo niot re- 
lations for a two-fluid shock (see also I Vinkl 120081 ). 

Here we explore the Rankine-Hugoniot relations for 
a two-fluid shock further. We show that for a given 
fractional cosmic-ray pressure w(= P C i/Ptat)> there 
is a unique cosmic-ray energy escape flux, e esc (= 
F cr / t^pqVs), associated with it that only depends on the 
overall Mach number of the shock. 

In the next section we present the derivation of the re- 
lation between w, e csc and post-shock plasma tempera- 
ture. In Sect.[3]we discuss this relation and its limitations 
in the context of two-fluid models, and cosmic-ray accel- 
eration models, and we use the relations to derive the 
cosmic-ray pressure content for the northeastern region 
of RCW 86 and for the newly measured plasma tempera- 
ture of the young L arge Magellanic Cloud SNR 0509-67.5 
(|Helder et alJl2010h . 



2. THE RELATION BETWEEN ESCAPE AND PRESSURE 
OF COSMIC RAYS 

Efficient particle acceleration by shock fronts leads 
to a shock structure that deviates sig nificantly from 
a normal "one fluid" shock fe. g. iDrurv fe Voelkl 



a no r mal one fluid s h ock I e. g. IDrurv &; Voelk 
19811: IBerezhko fe Ellisonl fl999t IBlasi et alJ 1200a 



Vladimirov et alJ 12 008: Dr urv et ; al.l l2009t IKang et al l 
2009t IReville et all 120091 ): the particles diffusing ahead 
(=upstream) of the shock, form a shock precursor that 
pre-compresses and slows down the gas flowing into the 
shock. The pre-compression caused by the precursor 
adiabatically heats the gas. The Mach number at the 
shock is, therefore, reduced with respect to the overall 
Mach number, as the shock velocity is reduced and the 
gas pressure upstream of the shock is increased with 
respect to a shock without a precursor. Additional 
heating of the gas may occur in the precursor due to 
non-adiabatic processes such as Alfvenic h eating (e.g. 
iVladimirov et all 120081: iCaprioli et all 120081 ). In such a 
multi-fluid system, the shock that heats the plasma is 
called the subshock. In the limit of a one fluid gas the 
subshock is identical to the shock. 

Her e we follow a differe n t app ro ach then | Blasi et all 
(2005 ): IVladimirov et all (120081) : IKang et all (I2009D : 
IReville et all (|2009f ) in that we treat the whole sys- 
tem only thermodynamically, using a two-fluid approach, 
with the two components representing the thermal gas 
and a gas of accelerated particles (cosmic rays). For 
the moment we neglect the possible influence of non- 
adiabatic heating in the shock precursor due to inter- 
actions between the gas and accelerated particles. Our 
appro ach is re miniscent of the work by IDrurv fe Voelkl 
(119811) (see alsolAchterberg et al.lll98llVoelk et al.1119841: 
IDrurv et al.ll2009D and is an extension of the work pre- 
sented bv lVinkl (|2008l ) and IHelder etHI (|200l . 



2.1. The Dependence of Post- shock cosmic-Ray Pressure 
on the Overall and Subshock Compression Ratio 

Our starting point is the basic Equations expressing 
conservation of mass, momentum and energy flux. We 
evaluate these expressions for three distinct regions: (0) 
far upstream, where the presence of accelerated particles 
can be neglected, (1) in the shock precursor, just up- 
stream of the subshock, and (2) behind (downstream of) 
the shock. Mass flux conservation gives: 



(1) 



with v the velocity of the gas with respect to the shock. 
Note that vq — V s , the shock velocity in the observer's 
frame of reference. Momentum flux conservation can be 
expressed as: 



Pq + poK 2 = Pi + pivi = p 2 



P2V2, 



(2) 



with P the pressure. 

At this point it is convenient to introduce the Mach 
number Mq for the shock structure far upstream, 



1 PoVs 2 
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with 7 g the adiabatic index of the thermal particles. In 
addition we introduce the compression ratios across the 
different regions 0,1, and 2: 
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We now use the assumption that the thermal pressure 
in the precursor is only due to adiabatic heating i n the 
precursor, i.e. P l th = P oX 1 s ( c .f. IDrurv et~aHl2b 091 and 
that across the subshock the pressure ass ociated with the 
accelerated particles do es not change ( see IDrurv fc Voelkl 
1981; Achterbcrg 2004). The cosmic-ray terms therefore 
cancel each other in Equation ([2]), when considering the 
momentum flux from region 1 to 2. In other words the 
non-thermal particle pressure at the subshock is only rel- 
evant for reducing the Mach number, but leads otherwise 
again to the standard, one fluid, relation for shock com- 
pression with a reduced Mach number given by: 
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Combining Equations ©, (O, and ((H) one finds for the 
total downstream pressure P2 and the thermal pressure 
Pth- 
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with the subscript cr referring to the cosmic rays. We 
have introduced here the symbol w describing the frac- 
tion of t he downstream pressure contr ibuted by cosmic 
rays ( c.f . Ivlnl 12001: IHelder et alJl2009h : 
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Post-shock temperature, cosmic-ray pressure and cosmic-ray escape. 
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Fig. 1. — Left: The relation between energy flux escape (e e sc, Equation. I15j l and downstream fractional cosmic-ray pressure (w, Eq. [9jl 
for an effective adiabatic index for the cosmic rays of -y cr = 4/3, and for various values of the upstream Mach number (Mo). Right: the 
relation between overall compression ratio (xi2, solid lines) and subshock compression ratio (x2, dashed lines). 

the precursor and the downstream fractional cosmic-ray 
pressure for a given upstream Mach number Mq. 

2.2. The Escaping Energy Flux 

In order to determine the escaping energy flux carried 
away by particles diffusing away far upstream, we need to 
use the expression for conservation of energy flux across 
the shock, but with a modification in order to express 
the fact that energy flux may be lost from the system 
(e.g. iBerezhko fc Ellison! ll9M ): 
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w = P cr/ P tot Ql 

Fig. 2. — Left: The same as the left panel of Figure QJleft), but 
now for an adiabatic index for the cosmic rays of 7 cr = 5/3. 

Using p\v\ = poVg /xi we can now derive an expres- 
sion for the fractional cosmic-ray pressure for a given 
Mach number Mq\ 



0) 



with u — P/(7 — 1) the internal energy, and 

r.r 



,)\p^}v s , 

(12) 



(13) 



(1-X?)+7 9 M 2 (1- i 



the escaping cosmic-ray energy flux, normalized to the 
total kinetic energy flux of the shock. Note that the 
escaping energy flux can only be taken out of the kinetic 
energy flux, as this is the only source of free energy. If 
we would have considered radiative losses, then the factor 
(1 — e C sc) should have been in front of Pq as well, as the 
usptream th ermal energ y can also be radiated away. 

Following |VE3 (p008h : iHelder et all (pOOl ) we intro- 
duce for convenience 



For high Mach number shocks this Equation simplifies 
to: 

I _ J_ _ 

W = t— = — . (10) 
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The compression ratio of the subshock is given by the 
standard shock relation: 



7s - 1 7cr - 1 7s-l 

Reordering the terms and using Equation (JB]) gives 
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(15) 
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with All being given by Equation 

Equations ([S]),® and (jlll) together show that there is 
a one to one relation between the compression ratio in 



This Equation completes the thermodynamic relation 
between the precursor compression ratio xi an d down- 
stream non-thermal pressure and overall energy flux es- 
cape. The resulting relation between e csc and w can be 
seen in Figure [TJ For high Mach number shocks the sec- 
ond and third terms can be omitted. 
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In the limit of 7 cr = 5/3 and Mo, Mi — > oo one can 
show that the relation between e eS c and w is well ap- 
proximated by: 



£csc _ A 1_\ 

W \ X\2> 



(16) 



Since realistically 7 cr < 5/3 (see Discussion), Equa- 
tion (|16l) serves as an upper bound on the escape flux. 

For low Mach number shocks there is a critical Mach 
number, M ~ 6 below which cosmic-ray acceleration 
cannot be thermodynamically supported for 7 cr =4/3, 
as can be seen in Figure [TJ For 7 cr = 5/3 (Figure [5]) 
cosmic-ray acceleration can be supported for lower Mach 
numbers, but there is still a limit, since e CS c < for 
M < 2.5. 

To summarize: taking xi as the principal variable, and 
Mo as an input parameter, one can calculate the sub- 
shock Mach number Mi (Equation from which the 
subshock compression ratio follows (Equation (fTTj) . and 
therewith the overall compression ratio X12 = XiX2- This 
can then be used to calculate w (Equation ([9])) and e eS c 
(Equation ([T5])'). We illustrate this in Figure [TJ 

2.3. The Maximum Overall Compression Ratio 

Note that although for a given w the escape flux e csc 
can be calculated, the opposite is not true. The reason 
is that e eS c has a maximum. In fact, mathematically 
the curves beyond the maximum declines further then 
indicated in Figure [TJ (left), but this would correspond 
to an unphysical compression factor of xi < 1 at the 
subshock. The peak value in the overall compression 
ratio corresponds to a maximum in e CS c • The peak value 
can be found by differentiating the expression for the 
total compression ratio (see Equation (flTTO 
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with respect to xu and setting dxxijdxx = 0. This gives: 
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Inserting this in Equation (lll[) using Equation ([5]) shows 
that X12 reaches a maximum for a subshock compression 
ratio of 



X2,( 



7s 



= 5/2, 



(19) 



7s -1 

with the numerical value valid for 7 9 = 5/3. 

2.4. The Downstream Temperature 

From Equations © and © or (|7J it is possible to find 
an expression for the downstream temperature: 

k B T 2 =P2,th/n 2 = (20) 
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with k B the Boltzmann constant and p the mean mass 
per particle in units of the proton mass m p . This expres- 
sion should be compared to the temperature expected 



behind a strong single-fluid shock: 
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(21) 



which can be obtained from Equation ([2T3]) by setting 
X12 = 4 (valid for a strong shock with j g = 5/3), w = 
and Mo — > oo. 
It is useful to define the ratio, 



k B T2 



(22) 



between the downstream temperature in the presence of 
cosmic-ray acceleration and the expected temperature for 
a pure gas shock, as this is a q uantity that can be relate d 
to existing measurements (c.f. Holder et al. 2009, 2010). 
The behavior of /3 is shown in Figure [3] (left). Note that 
Equation (|2"0"|) . like the expression for w (Equation , 
only relies on momentum conservation. As a result there 
is a unique relation between /3 and w. Since w does not 
depend on the adiabatic index 7 cr of the accelerated par- 
ticles, P does not depend on 7 cr , and, therefore, /3 does 
not depend on the energy distribution of the accelerated 
particles. 

In addition, we define the ratio a — k B T2/k B To be- 
tween downstrea m and upstream t e mperature, a s it is a 
quantity used bv lBlasi et all (|2005D ; IDrurv et all (|2009D . 
The expression for this quantity is: 



k B T p V? X12 



(23) 



(1-w)— h +7g M 2 (l- — 



1 



X1X12 



Xl2 

7s+l 
Xl 



Xl2 

7sM 2 (l-^ 

V X2 



The behavior of a as a function of w is shown in right- 
hand panel of Figure |3l In addition we show in Figure [4] 
the downstream Mach number Ma. This shows that for 
high w the dowstream Mach number rapidly approaches 
M% = 1. For X2 = 1, i-e. a continuous shock, we have 
M2 = 1, which leads to an unstable situa tion, as was 
already established by iDrurv fc Falld 



2.5. The Potential Influence of Magnetic Field 
Amplification 

As mentioned in the introduction, there is now ob- 
servational evidence for magnetic field amplification by 
shocks in young SNRs. In our two-fluid approach we ig- 
nore the effects of magnetic fields on the flow parameters, 
although we hope to come to this issue in the near fu- 
ture. For now we offer an assessment of the magnitude 
of magnetic amplification on the cosmic-ray dominated 
shocks. 

We start by noting that the downstream mag- 
netic field in several y oung SNRs is consistent wit h 
(B^/8n)/p V s 2 « 1 % (|Voelk et al.1 [20051: fVlnkl [2001 . 
The question now is, what is the perpendicular magnetic 
field pressure on the subshock compared to the ram pres- 
sure, i.e. what is (B 2 ± /8tt)/(pivD? This can be calcu- 
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Fig. 3. — Left: The ratio /3 between downstream (post-shock) temperature in a cosmic-ray dominated shock and a normal, high Mach 
number gas shock as a function of fraction downstream cosmic-ray pressure w. Right: The ratio <r between upstream and downstream 
temperature as function of w. 

supernova, magnetic fields as high as 50 G have bee n in- 
ferred, and B 2 /(8tt) ~ 0.1P cr (e.g. lTatischel[2009l) . At 
such a level the magnetic fields may start to become dy- 
namically important. 

Apart from amplifying the magnetic field upstream of 
the shock, cosmic rays in the precur sor may also give 
rise to non-adiabatic gas heating (e.g. IVladimiro v et al.l 
l2008t iCaprioli et al.ll2009t) . The effects of non-adiabatic 
heating will be similar to having a lower Mach number 
shock, i.e. non-adiabatic heating in the precursor limits 
the maximum overall compression ratio and limits the 
maximum possible fractional cosmic-ray pressure in the 
downstream region. 




0.4 0.6 

w = P =r/ P total 

Fig. 4. — The downstream Mach number M2 as a function of 
the fractional cosmic ray pressure w. 



2.6. The Entropy Change 

The en tropy jump across a shock is given by the rela- 
tion (e.g. IZel'dovich fc Raized Il966h : 
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Since B\ ± < B 2 we can say that observations of young 
SNRs indicate that 
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From Figure [T] (right) we can see that for a large range 
of w the factor X12/X2 1S smaller than one. Around the 
maximum of X12 we have X2 = 2.5. So in order to have 
(B%/8n)/p V? > 10% one needs xu £ 158, which is 
larger than for any of the models depicted in Figure [TJ 

However, there has been a debate whether B 2 oc poV 2 
(IVoelk et alJl2005h . or B 2 oc p V* (|Vinldl2008t iTatischeffl 
2009). For the young SNRs the precise proportionality 
does not change our conclusion very much, but in very 
young SNRs and/or for radio supernovae this may be rel- 
evant. In particular for SN 1993J, which was bright radio 



This neglects the en tropy increase assoc iated with the 
accelerated particles. I Brown et al.l (|1995[ ) calculated the 
entropy of non-thermal distributions, which, not surpris- 
ingly, always have entropy values below a thermal dis- 
tribution. Moreover, they normalize the distribution to 
the number of particles involved. However, in reality the 
cosmic-ray component may contribute a dominant frac- 
tion of the pressure, but the number of particles is always 
much smaller than the number of thermal particles. So 
neglecting the entropy of the cosmic rays is a good ap- 
proximation. 

Figure [5] (left) shows that, when the fractional cosmic- 
ray pressure w increases, the jump in entropy decreases. 
The right hand panel is more interesting as it shows that 
for certain values of the escaping energy flux e esc there are 
two possibilities for the jump in entropy. The maximum 
value of e eS c as a function of entropy change AS again 
occurs for the value for which X12 has a maximum. One 
can speculate that whenever there are possible solutions 
to obtain a certain value for e GSC the nature chooses the 
one that offers the highest entropy jump. On the other 
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Fig. 5. — The increase in entropy from far upstream to downstream. Left: the entropy increase as a function of downstream fractional 
cosmic-ray pressure w. Right: idem, but as a function of escaping energy flux £ MC . 



hand the escape has to be facilitated by the details of 
the acceleration process and particle spectrum, so one 
should not discount the low entropy branch too easily. 

3. DISCUSSION 

The solutions that we presented here for the two-fluid 
approach to shock with a cosmic-ray component give us 
a handle on estimating the effects of cosmic-ray accel- 
eration on the post-shock plasma temperature, and in 
addition allow us to estimate the cosmic-ray escape. It 
may not be too surprising that there is a thermodynamic 
relation between fractional cosmic-ray pressure w and es- 
cape. The reason is that in order to have a shock jump 
one needs either viscous shock, or one needs to lose en- 
ergy from the system. The problem with a cosmic-ray 
component is that the accelerated particles do not cause 
a jump in entropy. Therefore, the increased overall shock 
compression has to be facilitated by energy escape. 

This situation will change once additional heating in 
the precursor caused by cosmic-ray heating is intro- 
duced. This heating is presumably the result of dissipa- 
tion of large-amplitude magnetohydrodynamic waves in 
the cosmic-ray precursor (Alfvenic hea ting) ; these waves 
are the result of cosmic-ra y streaming dVladimirov et all 
l2008tlCaprioli et al.ll2009[) . The most important effect of 
heating in the precursor will be that the Mach number 
at the subshock Mi will be decreased. So we expect that 
including precursor heating will resemble the solutions 
for lower Mach number shocks. However, the physical 
details of cosmic-ray heating are not well known, and an 
in depth study of cosmic-ray heating is beyond the scope 
of this paper. Note that in principal additional heating 
can be incorporated in the Equations, by parameterizing 
it in relation to the upstream Mach number and w. 

Another limitation of the two-fluid approach is that 
it assumes a steady state situation and a plane parallel 
geometry. In a non-steady state situation the shock rela- 
tions may be influe nced by shock acceleration in the past 
(jDrurv et alJ IT995). Energy flux conservation (Equa- 
tionll2p may therefore be violated. Moreover, the highest 
energy particles may take such a long time between two 
consecutive shock crossings that the shock velocity has 
appreciably slowed down. An early discussion of time 
dependent effects, as well as the influence of spherical 



expansion, can be fou nd in iDrurv et al.l (|1995l) (see also 
iBerezhko et alJ 119941). Recently non-s teady situations 
were investigat ed by Kang et all (120091) . using a kinetic 
approach, and iSchure et al.1 (|2010l ). using Monte-Carlo 
simulations of test particles coupled to hydrodynamic 
simulations of spherically expanding SNRs. 

In non-steady situations high compression ratios can 
be reached, even in the absence of escape of cosmic rays 
from the system, because in the precursor there is a larger 
energy flux associated with particles diffusing in the up- 
stream direction than in the downstream direction. In 
a steady state situation, this asymmetry is caused by 
escape of the highest energy particles. In a non-steady 
state situation, in which the cosmic- ray particle popula- 
tion is building up (Ka ng et"aLll2009f ) , this asymmetry is 
caused by the fact that the flow in the upstream direction 
is from a more energetic particle population that left the 
sub-shock more recently. It would be interesting to in- 
vestigate whether one can model non-steady state shock 
acceleration in the two fluid approach by broadening the 
definition of e esc so that it includes this flux asymmetry. 

Note that pressure equilibrium (Equation [2]) is proba- 
bly a good approximation even in non-steady state situ- 
ations. As a result, the relation between k^T and w is 
expected to be valid even in non-steady state situations. 

3.1. Thermodynamics versus Acceleration Properties 

In the present model the quantities \2, X12, w, P, and 
a are independent on 7 cr , and encompass the solutions 
for particle accelerat i on based on nonlinear kinetic mod- 
els (IBlasi et alJl20CI5t iVladimirov et al.ll2008t iKang et aD 
200a iReville et alJl2009j ). in which the same Equations 
of the mass flux and the momentum flux conservation are 
adopted. Typically the kinetic models depend on two free 
parameters: the injection efficiency and the maximum 
momentum, beyond which particles will escape from the 
system. 

As can be s een in Fig. El our formalism agrees with 
the results of IBlasi et al.l (|2005f ). but the solutions by 
IBlasi et all (|2005| ) all cluster around w m 0.9 and X12 
around the maximum possible compression factor. The 
two-fluid approach, however, allows for a broader range 
of solutions. 
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Fig. 6. — Comparison of the two-fluid model with Mo = 100 with the kinetic model of Blasi ct al. (2005) (crosses). The kinetic model 
was used with p m ax = 10 J mc, § = 3.5, and for 20 logarithmically spaced velocities in the range of 10 — 10 5 kms -1 . Left: cosmic-ray escape 
version post-shock partial cosmic-ray pressure. Right: overall and subshock compression ratio. 



In the two- fluid model e esc depends on 7 cr , which is 
essentially a free parameter. In the kinetic models e CS c is 
determined self-consistently by the overall hydrodynamic 
structure of the shock, including the dynamical effects of 
the accelerated particles. The energy distribution of the 
accelerated particles determines then 7 cr = 1 + P CT /u CI , 
where P CI is the cosmic-ray pressure and u CT the cosmic- 
ray energy density. Hence the kinematic models predict 
that the e esc — w diagram should deviate from ours with 
7cr = 4/3, and should lie between our curves of 7 cr = 4/3 
and 5/3. In practice, however, the cosmic-ray adiabatic 
index as obtained by the kinetic models is very close to 
7cr = 4/3 for p m ax = 10 5 mc, as can be seen in Fig. [51 
For p m ax = 10 2 mc we found that a modest increase to 
7 cr = 1.4 resulted in a good match between the two-fl uid 
solutions and the kinetic model of iBlasi et al.l (|2005[ ). 

Interestingly, the two-fluid solutions presented here in- 
dicate that a higher energy flux escape is necessary for 
7 cr > 4/3, but from the point of view of the spectral en- 
ergy distribution the energy escape flux is more difficult 
to achieve for spectra with spectral energy indices T > 2 
{N{E) cx E~ r ), corresponding to 7 cr > 4/3. The reason 
is that for T > 2 most of the cosmic-ray energy is con- 
tained by mildly relativistic particles, whereas most of 
the escaping particles will be near the maximum of the 
energy distribution (p m ax in the Blasi model). Hence, 
the escaping particles carry away only a small fraction 
of the internal energy. For T < 2 most energy is indeed 
concentrated around p m ax and the necessary escape flux 
is easily generated. 

So from a thermodynamic point of view, a high frac- 
tional cosmic-ray pressure w, requires large values of e e so 
which in turn is best achieved with 7 cr « 4/3. 

3.2. Application to Observations 

One of the motivations for this work has been the 
measurements of the downstream te mperature of SNRs 
with reasonably well known velocities (fHelder et al.ll2009l 
2010). These temperatures were measured using the 
broad component of the Ha line. This component is 
due to charge exchange between neutral hydrogen atoms 
entering the shock and the downstream population of 



shock heated protons. Hence, the width of the line is 
caused by thermal Doppler broadening and reflects the 
downstream proton temperature. There is some ambi- 
guity as to how to relate the proton temperature to the 
overall downstream plasma temperature, as the different 
plasma constituents (electrons, protons , helium, other 
ions) may not be in thermal equilibrium. However, the 
proton temperature is always expected to be within a 
factor of 2 of the mean plasma temperature. 

Using Equation (f20|) (see also Fig. [3]) one can easily 
estimate from a measured temperature and shock veloc- 
ity what the downstream fractional cosmic-ray pressure 
is and what the required escape flux is. The only ambi- 
guity that is left for estimating e CS c is what the effective 
cosmic-ray adiabatic index is. However, 7 cr is not impor- 
tant for determining w. For the SNRs under considera- 
tion the shock velocities are well in excess of 1000 kms . 
For a typical sound speed in the interstellar medium of 
10 kms -1 , we have Mo > 100. This means that the high 
Mach number approximation is valid as long asw< 0.9. 
This appears to be the case for the SNRs considered be- 
low. 

For the northeastern r egion of the TeV g amma-ray 
emitting remnant RCW 86 Held er et al.1 ([20091 ) measured 
a downstream temperature of k^T p = 2.3 ± 0.3 keV for a 
measured shock velocity of V s — 6000 ± 2800 kms -1 . 
Nominally this corresponds to (3 — 0.055 (see Equa- 
tion ([20]l). but given the systematic uncertainties and 
the ambiguity due to non-equilibration of temperatures j3 
could be as high /3 = 0.31. Using now Equations (j9|) and 
(|15|) in the limit for high Mach numbers the measured 
values of j3 correspond to downstream fractional cosmic- 
ray pressures and escape fractions of w = 0.81, e osc = 
0.59 and w = 0.51, e osc = 0.20, for 7cr = 4/3 and 
/3 = 0.055 and /3 = 0.31 respectively. For -f cr = 5/3 
this is e csc = 0.72 and e csc = 0.38, respectively, with w 
unchanged. 

F or the young Large Magellanic Cloud remnant 0509- 
67.5 lHelder et al.l (|2010l ) determined the post-shock tem- 
peratures in two regions. The most constraining mea- 
surement was for the southwestern region for which 
ksT p — 28.7 keV, for a conservative velocity estimate 
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of V s — 5000 kms^ 1 , corresponding to ft = 0.58. This 
translates into a downstream fractional cosmic-ray pres- 
sure of w — 0.29, e osc = 0.06 for 7 cr = 4/3 or e csc = 0.19 
for 7„. = 5/3 . 

Interestingly, these new estimates of w based on Equa- 
tio ns (|9| and (|15l) are not fa r from the lower limits given 
bv lHelder et alTpOOl l20l"oh . The reason is that the less 
constraining relations used by iHelder et al.l ((20091 12010ft 
allow in principle for a higher fractional cosmic-ray pres- 
sure by reducing the cosmic-ray escape flux. The rela- 
tions derived here do not allow for this possibility. More- 
over the relation that was used to determine the lower 
limit by IHelder et all (|2009l 120101 ) is a limiting case of 
the shock relations presented here (see Equation (fT6l) ). 



3.3. Do Very Efficiently Accelerating Shocks Exist? 

This brings us back to one of the principle issues cur- 
rently discussed in cosmic-ray physics, namely do very 
efficiently accelerating s hock exist? The pos t -shoc k tem- 
peratures measured by IHelder et al.l ([200 9. 2010) indi- 
cate high values for the cosmic-ray pressure, but not as 
high as usually found by particle acceleration models, 
where the efficiency is often found to be close to 90%, 
with total compr ession factors as high as X ]2 = 70 ( e -%- 
iBlasi et aD 120051 ). It has been argued by IDrurv et all 
<|2009f ) that such extreme conditions may exist in the 
Tominent TeV gamma -ray source SNR RX J1713. 7-3946 
Aharonian et al.ll2004l ). Apart from high gamma-ray lu- 
minosity this idea is based on the lac k of detectab l e ther - 
mal X-ray emission. According to IDrurv et al.l ()2009[ ) 
this could be due to a very low post-shock plasma tem- 
perature, possibly as low as low as six times the up- 
stream. Alternatively, RX J1713. 7-3946 may evolve in 
a low density medium, resulting in a low thermal lumi- 
nosity, as the thermal luminosity scales with the square 
of the density. The Equations derived here allow for 
ve ry low downstream temperatures. However, the claim 
bv IDrurv et~aT1 (|2009| ) that the downst ream temperature 
may only be six times the upstream temperature seems 
unrealistic, as Equation ([9]) and (|15j) show that this is 
only true for unrealistically low Mach numbers. For more 
realistic Mach numbers for young SNRs, say M > 75, 
the downstream temperature is at least 80 times the up- 
stream temperature. 

It could be that the young SNRs that we observe are 
already past their prime as cosmic-ray accelerators, al- 
though the cosmic-ray escape in RCW 86 appears still 
high enough to explain the cosmic-ray energy production 
of SNRs. Specifically it has been argued that in the very 
early SNR phase/supernova p hase cosmic-ray accelera- 
tion may be very efficient (e.g. iBell fc Lu cek 2001; iVinkl 



2008) , see for example the case of SN 1993J (iTatischefl 

2009) . Note that a steady escape of cosmic rays dur- 
ing the evolution of a SNR is favored over a more sud- 
den release of cosmic-rays at the end of the SNR evolu- 
tion, as in the latter case adi abatic losses have lowered 
the energy of the cosmic-rays (jVoelk et alJll984j ). If one 
wants to explain the cosmic-ray energy density in the 
Galaxy by cosmic-ray escaping from young SNR shocks 
(i.e. e csc = 5%, see introduction) then on average over 
the lifetime of the SNR a fractional cosmic-ray pressure 
is necessary of at least w ~ 30% for 7 cr =4/3. 

4. SUMMARY AND CONCLUSIONS 

We have shown that for a two-fluid steady state shock 
consisting of a thermal gas plus a cosmic-ray component 
there exist a unique relation between the downstream 
fraction of pressure provide by cosmic rays (w) on the 
one hand, and the energy escape flux normalized to the 
incoming free energy flux (e csc ) on the other hand. This 
relation depends on the assumed effective adiabatic in- 
dex of the cosmic rays and on the overall Mach number. 
For Mach numbers Mo < 6 no cosmic-ray pressure is al- 
lowed. This could be of interest for particle acceleration 
in clusters of galaxies, since internal cluster shocks tend 
to be low Mach number shocks, whereas the accretion 
shocks or cluster formation shock are high Mach number 
shocks. 

This relation can be used to determine w and e esc based 
on measurements of the downstream temperature of the 
plasma for a known shock velocity. Using this relation we 
show that the lower limits on w inferred for shock regions 
in the supernov a remnants RCW 86 and SNR 0509-675 
bv lHelder etUl (pOOl l2010h are in fact very close to the 
actual value for w, but we can now determine also e csc , 
which is at least e esc = 0.2 for the northeastern region of 
RCW 86, but only 0.067 for SNR 0509-675 (but it may 
be 0.19 in the unlikely case that the cosmic rays have an 
effective adiabtaic index of 7 cr = 5/3). 

In order to explain the cosmic-ray energy density in 
the Galaxy by cosmic-rays escaping directly from SNR 
shocks, instead of released at the end of the SNR life, one 
needs an escape flux of e C sc ~ 5%, which requires that on 
average 30% of the downstream SNR pressure should be 
supplied by cosmic rays. 
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